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CO ■ We investigate the possibility of the identification of TeV physics models including 

vO ' WIMP dark matter at the International Linear Collider. Many TeV physics models 

contain a WIMP dark matter (x°) and charged new particle (x ) which interacts with 
r — . the WIMP dark matter via the vertex x ± X°W zf: . Through Monte Carlo simulations, 

we study the process, e + e~ — > x + X~ ~~* X°X°W + W~ , because the signal contains the 
^J , fruitful information of the model. We show that, in particular, the distribution of the 

X production angle is the powerful probe in the TeV physics model search. 



1 Introduction 



About 23 % energy density of the present Universe is made up of dark matter and it plays a 
key role in the large structure formation. However, we don't know what it is. Dark matter 
candidates don't exist in the standard model (SM). 

Weakly Interacting Massive Particle (WIMP) is one of the most plausible candidate for 
dark matter. It is neutral, massive, and sufficiently stable particle and such particle naturally 
explains the observed dark matter abundance, SloMh 2 ~ 0.11 [I]. Many models of physics 
beyond the SM including WIMP dark matter have been proposed. 

The TeV physics could be produced at the Large Hadron Collider (LHC) and the Inter- 
national Linear Collider (ILC). The LHC is now operating, but the precision measurement 
of events including dark matter in final state is difficult at the experiments. On the other 
hand, it would be possible at the ILC. The measurements is important for not only dark 
matter physics but also identification of the TeV new physics model. 
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Models 


Particles 


Spins 


Inert Higgs like model 


(xS, xV 


(0,0) 


SUSY like model 


(xi, x° F ) 


(1/2, 1/2) 


Little Higgs like model 


(Xv, X"v) 


(1,1) 


(1,0) model 


(Xy, X°s) 


(1,0) 


(0, 1) model 


(x^,x u v0 


(0,1) 



Table 1: Classification of spin configuration of new particles. 



In this study, we investigate the possibility of general WIMP search and new physics 
model identification studying e + e~ — > x + X~ ~~> X°X°W + W~ process at the ILC, where x° 
and x^ denote the WIMP and the charged new particle, respectively. Using the process, 
we have already studied about the possibility of precision measurements of the Little Higgs 
model [2]. We have shown that we can extract fruitful information from the signal of 
the process: The masses of x^ an d X° can be determined from the edges of the energy 
distribution of the reconstructed W bosons. It is also possible to confirm that the spin of 
X and the vertex structure using the angular distribution of the x P a i r procudion and 
the polarization of W^. The gauge charges of the x^ boson could be measured using a 
polarized electron beam. 

It is worth noting that such process exists in many TeV physics models which explain 
the smallness of electroweak scale , because WIMP candidates would interact with the SM 
particles weakly Thus, the process plays a key role of the general WIMP search. 



2 Benchmark models and representative points 

We concentrate on the WIMP dark matter (x°) and the charged new particle (x^) which 
interacts with dark matter. The interaction vertex, x X°W^, is also exist in many TeV 
physics models which contain WIMP scenario. All spin combinations which has the vertex 
is written in Table [TJ keeping charge, gauge, discrete symmetry and Lorentz invariance. 

In this study, we investigate the Inert Higgs doublet like model (IH) [3] , supersymmetric 
like model (SUSY) |3], and Littlest Higgs like model (LHT) [5]- [TTJ as benchmark models of 
X^ spin 0, 1/2, and 1, respectively . The crucial difference from the (1,0) and (0,1) models 
in Table Q] only appears in the what relates with the x ± X°W /=F vertex (e.g. the shape of the 
energy distribution of W bosons). 

At representative points, we take the same new particle masses and production cross 
section in each models in order to investigate the separation possibility of the TeV new 
physics model using the same masses and number of events. Actually, when we discuss the 
features of new particles which observed at future colliders, we compare the signal with the 
sample models in the same way. The parameters at the representative points are summarized 
in Tabled 

Since the mass difference of new particles is determined in littlest Higgs model explicitly, 
the mass differences at representative points are adjusted to coincidence with Littlest Higgs 
model. We take the coefficients of vertices to consistent with a parameter in each models 
which realizes the mass difference. However, in order to adjust the amplitude of cross section 
to each other, the over all of neutral couplings (e.g. x + X~ %) are normalized. 
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m x ± [GcV] 


m x o [GcV] 


Cross section 


Br(x ± -4 X^W*) 


y/a = 500 GcV 


231.57 


44.03 


40, 200 [fir 1 ] 


~ 100% 


V^ = 1 TcV 


368 


81.9 


40, 200 [ftr 1 ] 


~ 100% 



Tabic 2: Representative points in this study 



Detector 



Performance 



Coverage 



Vertex detector 
Central drift chamber 
EM calorimeter 
Hadron calorimeter 



5 b < 5 9 lO/p/3 sin 3/2 6 H 
WPt < 5 x 10~ 5 (GeV/c)- 1 

a B /E = l7%/VE®i% 

er B /£ = 45%/V£©2% 



cos6»| < 0.93 
cos 6»| < 0.98 
cos(9| <0.99 
costfl <0.99 



Tabic 3: Detector parameters used in our simulation study. 



In this paper, we show the study at y/s = ITeV with an integrated luminosity of 500 
fb -1 . The study at y/s = 500GeV is written in [12]. 



3 Simulation tools 

In our study, both signal and background events have been generated by Physsim JT3] . The 
initial-state radiation and beamstrahlung have been included in the event generations. We 
have ignored the finite crossing angle between the electron and positron beams. In the event 
generations, helicity amplitudes were calculated using the HELAS library [T3], which allows 
us to deal with the effect of gauge boson polarizations properly. Phase space integration 
and the generation of parton 4-momenta have been performed by BASES/SPRING [T5] . 
Parton showering and hadronization have been carried out by using PYTHIA6.4 [IB], where 
final-state tau leptons are decayed by TAUOLA [17] in order to handle their polarizations 
correctly. 

The generated Monte Carlo events have been passed to a detector simulator called JS- 
FQuickSimulator, which implements the GLD geometry and other detector-performance 
related parameters [18] . In the detector simulator, hits by charged particles at the vertex 
detector and track parameters at the central tracker are smeared according to their position 
resolutions, taking into account correlations due to off-diagonal elements in the error ma- 
trix. Since calorimeter signals are simulated in individual segments, a realistic simulation 
of cluster overlapping is possible. Track-cluster matching is performed for the hit clusters 
in the calorimeter in order to achieve the best energy flow measurements. The resultant 
detector performance in our simulation study is summarized in Table [3] 



4 Signal Selection 



->x°x 



We use events in which both W bosons decay hadronically, e + e~ — > x X 
X°X qqqQi as signal events. Thus, the main background processes are VF+VT 
etc. 



°w+w- - 

vvW+W~ 
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Physics model 


<T = 40fb 


a = 200 fb 


M x ± (GeV) 


IH 


364.3 ± 6.0 


366.4 ± 1.4 




SUSY 


370.6 ± 5.6 


368.0 ± 1.3 




LHT 


367.7 ± 4.0 


367.4 ± 0.9 


M x o (GeV) 


IH 


79.7 ± 4.8 


78.5 ± 1.2 




SUSY 


76.6 ± 6.5 


77.8 ± 1.3 




LHT 


78.0 ± 3.9 


78.7 ± 0.9 



Table 4: Measurement accuracy of masses of x ± and x°- 



In order to identify the two W bosons from x decays, two jet-pairs have been selected 
so as to minimize a x 2 function, 



X 



r 



IVKl 



tr M ff )>; v 



( rec M 



W2 



tr M^)>M w , 



(1) 



where I0C M W1 ( 2 ) 1S the invariant mass of the first (second) 2-jet system paired as a W 
candidate, tr Miy is the true W mass (80.4 GeV), and <7m w is the resolution for the W mass 
(4 GeV). We required x 2 < 26 to obtain well-reconstructed events. 

Since x° escape from detection resulting in a missing momentum, we have thus selected 
events with the missing transverse momentum mlss p T above 84 GeV. We have also selected 
events with a energy of W below 500 GeV. The number of remaining background events is 
much smaller than that of the signal after imposing all the cuts. 



5 Mass Determination 

The Masses of new particles, x° an d X , can be determined from the edges of the W energy 
distribution. After subtracting the backgrounds, the distribution has been fitted with a line 
shape determined by a high statistics signal sample. The fitted masses of x° and x^ are 
summarized in Table[4] The masses of x and X° W1U be determined with 0(1)% accuracy. 

6 Angular Distribution 

The production angle of x can be calculated with 2-fold ambiguity from the momenta of 
W bosons, assuming back-to-back production of x + and X- Figure [1] shows 2-dimensional 
histogram. To estimate the possibility of the new physics model identification, we compare 
the distribution with the " template" which is the distribution of high statistics signal samples 
using x 2 



2 analysis: 



^Data — -Vremp 
CSig+BG 



(2) 



where iVo a ta and iVremp are the number of the data and the number of the template at each 
bin, respectively. The <7si g +BG ^ s the error of the signal and backgrounds. The % 2 values 
are summarized in Table [5] It shows that the identification of these new physics models is 
possible by comparing of production angles of x ■ 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

(IH) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

(SUSY) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

(LH) 



Figure 1: The 2-dimensional histogram for the two solutions of the two solutions of the 
production angle at each models. 



a 


Model 


IH template 


SUSY template 


LHT template 


40 fb 


IH 

SUSY 
LHT 


1.2 

5.2 
3.7 


4.7 
1.2 
2.6 


3.6 
3.2 
1.2 


200 fb 


IH 

SUSY 
LHT 


1.0 
20.0 
13.7 


17.3 
1.0 
7.8 


11.0 
9.2 
1.0 



Table 5: Summary of x 2 /(the number of bins) values from 2-dimensional histograms of 
X production angle. 



the 



7 Summary 

We study the possibility of general WIMP search and new physics model identification using 
e + e~ — > x + X~ -^ X°X°W + W~ process at the y/s = 1 TeV ILC with integrated luminosity 
of 500 fb -1 . We have shown that the masses of new particles can be determined very 
accurately at the ILC in a model independent way. Furthermore, using x 2 analysis, the 
identification of new physics models is possible by comparing of the x production angles. 
Finally, we have also studied the threshold scan to separate the new physics model. It shows 
that the SUSY like case will be separate from other models. 
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